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ABSTRACT 


This report considers the motion of one spacecraft in a near-cir- 
cular orbit which is fastened to another with a flexible tether line. 

The equations of motion are obtained through a closed form solution to 
the approximate differential equation. A method is devised by which 
pulls can be made on the tether in a repeatable sequence so that the 
two spacecraft will remain in the vicinity of one another when there are 
differential drag forces acting on them. This method is then applied to 
the case of a Lunar Module/Apollo Telescope Mount (LM/AIM) tethered to 
an S-IVB Orbital Workshop (OWS). 
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A METHOD OF SOFT TETHER STATIONKEEPING 
SUMMARY 


This report considers the motion of one spacecraft in a near-cir- 
cular orbit which is fastened to another with a flexible tether line. 
The equations of motion are obtained through a closed form solution to 
the approximate differential equation, A method is devised by which 
pulls can be made on the tether in a repeatable sequence so that the 
two spacecraft will remain in the vicinity of one another when there 
are differential drag forces acting on them. This method is then 
applied to the case of a Lunar Module/Apollo Telescope Mount (LM/AIM) 
tethered to an S-IVB Orbital Workshop (OWS). 


I. INTRODUCTION 


With the advent of more sophisticated space missions, the need to 
maintain multiple spacecraft in close proximity to one another arises. 
Intuitively, it seems that there would be no problem caused by the 
extremely small differential forces acting on these spacecraft arising 
from small variations in venting thrusts, solar pressure and drag 
decelerations. However, as figure 1 shows, these forces, typical of 
those acting on two spacecraft of differing configuration and orienta- 
tion, can cause large separation distances after several days* time. 

For this particular figure, the difference in the drag deceleration 
between the two bodies was .312 x 10~®m/sec^. 

There is a variety of methods which may be employed to prevent 
the separation of spacecraft in orbit. They may be mechanically fixed 
together (hard docked) or they may be unattached and kept together by 
propulsive means. A compromise between these methods is to fix the 
bodies together with a tether line. Each of these methods has advan- 
tages and disadvantages depending on the goals of the particular mis- 
sion. For missions which are propellant-limited but require that two 
initially closely spaced vehicles have independent attitudes, the tether 
method looks attractive enough for a more detailed analysis. 

The tether method can be divided into two classes. One of these 
is characterized by a continuous tension in the tether. In the other 
class, the tether is taut at only discrete time points and is slack 
the rest of the time. This report analyzes the relative motion of two 



spacecraft for the latter class and presents a method for controlling 
their relative positions within desired bounds. No attempt is made to 
describe the methods which may be used for deployment of retrieval [1] 
of such tethered satellites. 

To make an analysis of the spacecraft motion, first, the differen- 
tial equation of motion of one spacecraft relative to the other is 
linearized and solved in closed form. This solution is then used to 
find a method by which the two spacecraft can be kept together. Finally, 
the physical characteristics of this method are described for the example 
of a Lunar Module/Apollo Telescope Mount tethered to an S-IVB Orbital 
Workshop. 


II. DERIVATION OF THE EQUATIONS OF MOTION 


Let one of the spacecraft be in a circular orbit of radius, R, and 
angular velocity, o), (figure 2). Set up a two-dimensional coordinate 
system in this spacecraft with the positive y-axis pointing outward 
along the radius vector and the x-axis in the opposite direction of the 
inertial velocity vector. Let the position vector in this coordinate 
system to a second spacecraft be r. If R is the radius vector from the 
center of the earth to the origin of the relative rotating coordinate 
system, then the differential equation of motion for the center of mass 
of the second body can be written: 


f = a - 2Zsx. f + g - cjx [cox(f + R) ] , 


where a is the sum of all the differential accelerations between the two 
bodies due to aerodynamic, propulsive and tether forces, -2cax f is the 
Coriolis acceleration, g is the gravitational acceleration, and the last 
term is the centrifugal acceleration. The last two terms can be com- 
bined in the following manner: 

LL(r + R) 

g = • , 

(r + R)2 

according to the inverse square gravitational law, where jj, is the gravi- 
tational constant for the earth. Since = M./R^> 

- _ -u)^(r + R) 

s • 

(r + R)^ 


2 


Since co is perpendicular to (r + R), 


and 


OK [ojix(rx5.)] = -w^(r + R) 


g - (jjx[co + (r +R)] = w^(r +R) 

\ (r 

The terra (r + R) can be approximated by R + y with an error of only 
about x^/R, where x and y are the components of f in figure 2. There- 
fore, 



g - ojx[oK(r + R)] = o)^(r + R) 


1 - 


(R + y)' 


= co^(r + R) I 1 - 


Z^3. 


(1+t) 


Expanding the second term in the brackets by the binomial expansion 
yields 


g - ojx[oK(r + R)] = oo^(r -h R) 




^ P (r + R) 


where ^ may be approximated by the unit vector in the y direction, j, 


R 


f = a - 2cax r + 3to^j, 


Separating this acceleration into x and y components yields 


x = ajj + 2(jjy 
y = a “ 2(jjk + 3w^y. 

y 
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For this report, it will be assumed that the difference in drag 
deceleration acts on the second body in the positive x-direction and is 
a constant, D. The Rand Corporation investigated the equations for a 
drag which is dependent on altitude [2]. However, this assumption is 
justified for very nearly circular orbits in a spherically symmetric 
atmosphere (with no diurnal bulge). During the time when there is no 
tension in the tether, these differential equations therefore become 


X = D + 2coy 

(1) 

y = -2ojx + Sco^y. 

(2) 


These equations may be solved analytically as shown in Appendix A to give 
X, y, X, and y as functions of time. 


where 


X = + c^t - Dt^ - 2b^ sin cot + 2aQ(l - cos cot) 


(3) 


2Co 2 d 

y = 7 — t - b^ cos cot + sin cot 

3to CO o o 


(4) 


X = c^ - 3Dt - 2 c±)q cos cot + 2coa^ sin cot 


(5) 


2D 


y + cob_ sin cot + coa^ cos cot, 

M o u » 


( 6 ) 


1 f. ^ 2d\ 

a = - ( y + — , 

O 0} V 0 GO ' 




c^ = 6coy_ - 3x , 
o *^0 o 


and Xq, yo, Xq, and yo are the conditions of position and velocity at 


time t = 0, 
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These equations are not restricted to use with tethered satellites. 
They may be used to describe the mption of any body in a near-circular 
orbit which remains fairly close (i,e,, much less than the radius of the 
orbit) to the origin of the rotating relative coordinate system. For 
example, they could be used for the coasts between burns of the terminal 
maneuvers of rendezvous and docking. Also, these equations are not 
restricted to describing the motion of only one body. By having a 
separate set of equations as (3) through (6) for each body, the motion 
of multiple spacecraft in low eccentricity orbits could be described. 

By subtracting their coordinates, the motion of all of the bodies could 
be described relative to one of them. This is a possible method of 
removing the restriction that one of the spacecraft be exactly in a 
circular orbit. 


III. SOLUTION TO THE STATIONKEEPING PROBLEM 


In the following discussion, the body located at the origin of the 
relative coordinate system will be called the reference body and the 
other spacecraft the tethered body. 

The equations of motion developed in the previous section cause 
the tethered spacecraft to follow a trajectory relative to the refer- 
ence body similar to that shown in figure 3. In this particular example 
the motion of the tethered body is initiated at the position x = 30 m, 
y = 4 ra, with a positive x and y such that its orbital energy is greater 
than that of the reference body. Superimposed onto an oscillation in 
the vertical and horizontal direction, the tethered body first experi- 
ences a drift to the right and then a drift to the left. The oscilla- 
tion in the vertical direction results as it moves between apogee and 
perigee. The oscillation in the horizontal direction is due to the fact 
that the perigee and apogee velocities are different. The drift to the 
right is caused by the semimajor axis of the conic of the tethered body 
being initially greater than that of the reference body and thus having 
a longer period. As drag decays the orbit; the semimajor axis decreases 
and the spacecraft then drifts to the left. The orbital energy lost 
through atmospheric drag must in some manner be regained, or the orbital 
period will continue to become less and less, thus increasing the separa 
tion rate. If energy is to be added to the tethered body by adding to 
the velocity by pulling on the tether, this must be done before it gets 
ahead of the reference body (since any pulls after it gets ahead will 
only decrease its velocity). It would also be highly desirable to be 
able to get on a transfer ellipse (such as indicated in figure 3 by the 
dashed line), and then when the spacecraft intersects the original 
trajectory, give a second impulse to reproduce the initial conditions. 

If this could be done, the process could be repeated over and over. 
Obviously, the correct impulses to do this could be found if the 
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magnitude and direction of the pull could be controlled. However, the 
direction is restricted to lie in the direction of the taut tether. 

Since, in general, this will not be the desired direction, some constraint 
must be placed on the initial trajectory which will insure that the 
direction of the tether be correct to initiate the transfer and obtain 
a repeatable sequence. In Appendix B, it is found that a sufficient con- 
dition to insure repeatability is for the motion of the tethered body to 
describe a trajectory which is symmetric with respect to the x-axis in 
the relative coordinate system. Figure 4 shows a typical symmetric tra- 
jectory. Each one of those showing symmetry can be characterized by 
two arbitrary parameters. Let these be the initial altitude difference, 
yi^, and the time, At^, spent from the initial point until a pull on the 
tether is required. For any combination of these two, there can be found 
an initial velocity which will cause the motion to be symmetric. As 
shown in equations (B21) and (B22) in Appendix B, the components of this 
velocity are independent of the initial tangential separation, x^^. How- 
ever, there is a restriction on x^^ to be positive (the tethered body 
behind the reference body). This restriction is made to make equation 
(B31) solvable with a positive AV^ since a negative AVi would be equi- 
valent to a push with the tether. 

In summary, the stationkeeping problem can be solved with a repeat- 
able sequence by placing the tethered body on a symmetric trajectory. 


IV. A PRACTICAL APPLICATION 


To obtain greater insight into the characteristics of the station- 
keeping problem and the use of symmetric trajectories as a solution, 
first a specific tether mission is considered. Then several of the 
variables pertinent to this mission are parameterized. Finally, the 
most attractive symmetric trajectory is chosen and described in more 
detail . 

An S-IVB orbital workshop (OWS) [3] space station is in a 260-. 
nautical-mile circular orbit. It is desired to keep a Lunar Module/ 
Apollo Telescope Mount (LM/ATM) [4] in the vicinity of the workshop by 
use of a tether line. As shown by the aerodynamic data in table I, the 
ATM will experience the greater drag deceleration due to orientation and 
ratio of drag to mass. The difference in drag acceleration, D, is about 
.312 X 10“® m/sec^. 

The variables which were considered as being important or charac- 
teristic of the trajectories are the initial velocity components, the 
transfer time, At 2 , the total cycle time, At^ + At 2 , and the amplitude 
of the perigee-apogee oscillation. The initial velocity components 
necessary to start the ATM on a symmetric trajectory can be obtained 
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from equations (B21) and (B22) in Appendix B. Since these components 
are functions only of the altitude difference at the beginning of the 
symmetric trajectory, y^i, and At^, they can be represented as contours 
on an abstract plane, y^^^, At^ (see figures 5 and 6). Similarly, the 
transfer time, Atg, characteristic for each S3mtmetric trajectory may 
be plotted on this plane for some specified value of x^^i, the initial 
tangential separation. In figures 7 through 10, this has been done 
where Xij^ » Multiple transfer times are available for any sym- 

metric trajectory due to the trigonometric part of equation (B33) . By 
making the same assumption on figures 11 through 14 show contours 

of equal total cycle time, At^ + At 2 . In this case, by multiplying by 
1/2 D, these become contours of equal AVi. For a cycle time of 4X, 

AVi is .003527 m/sec. To impart this velocity change to the LM/ATM, 
it is necessary to pull on the tether with a force of 1 pound for 
5.86 seconds. This impulse being typical, only a very light tether 
line is required— of the order of fishing line. The amplitude of the 
oscillation as the ATM moves between apogee and perigee may be derived 
as a function of y^j^ and At^ by using equation (B18). Contours of equal 
amplitude, ap will obey the equation 


11 


= ± a 


sin uAt, + — At,. 
2 10 


These contours are shown in figure 15. On an actual mission, At^'s near 
an integer number of orbital periods should be avoided because of large 
oscillations and because of velocities in the vertical direction. 

The family of solutions indicated in figure 15 as having zero ampli- 
tude in the vertical direction is particularly interesting. The solid 
line in figure 16 shows the appearance of one of these case. In these 
cases, both the workshop and ATM are in circular orbits,* the ATM orbit 
initially being slightly higher. Since the period of the ATM is thus 
longer, it will drift further behind the workshop. Drag gradually lowers 
the altitude of the ATM orbit, decreasing its period. After its altitude 
drops below that of the workshop, it begins to catch up. When it has 
come to within a certain distance, Xif> a tug is given which places it on 
a transfer trajectory to the initial conditions where a second pull 
recircularizes it. If x^i » (yii), then the transfer is very close to 
a Hohmann transfer making Atg about equal to half an orbital period. 


The appearance of the trajectory would be very similar if the two 
bodies were in elliptical orbits of equal eccentricity and coincident 
lines of apsides. 


7 



The Ati and the maximum distance between the ATM and workshop, 
which are obtained, and the Xi^, which is necessary, are plotted as 
functions of yii in figure 17. The necessary y^^ is always equal to 
-2 d/co. These cases are fairly insensitive to errors in both the timing 
and magnitude of the pulls on the tether. First, consider the pull 
which is made below the altitude of the workshop. If this pull is made 
somewhat early or late, the transfer back to the initial conditions will 
be performed slightly farther from or nearer to the workshop than the 
nominal. If the magnitude of the pull is not quite nominal, then the 
transfer will be to a somewhat different altitude. When the recircular- 
ization pull is made, the ATM will merely be on a different member of 
the family of zero oscillation cases indicated by the dashed lines in 
figure 16. Errors in the magnitude of the second pull as much as ± 10 
percent and in timing as much as from 3 to 5 minutes introduce only a 
small amount of oscillatory motion, as indicated in figures 18 through 
21 . 


V. CONCLUSIONS AND EXPANSIONS 


It is possible to conclude from the ground work set up in this 
report that it is feasible to. overcome the separation buildup between 
two satellites under the influence of differential drag forces by using 
a flexible tether line. This can be done by bringing the tether taut, 
imparting an impulse through it, and allowing the bodies to drift free 
before a second impulse is made. The sequence for pulling on the tether 
becomes repeatable if the motion of the tethered body appears symmetric 
relative to the reference body. The most practical of these symmetric 
cases seems to be a case which exhibited no oscillations. 

Several areas which need further work are being investigated. The 
reduction of the equations in Appendix B should be generalized to estab- 
lish the necessary conditions for a repeatable sequence.* Also, motion 
only in two dimensions has thus far been considered; this should be 
extended to motion out of the orbital plane. Errors in the timing and 
magnitude of the impulses should be further investigated, and some con- 
trol laws should be developed to maintain a nominal sequence. Before 
making any actual flights, particularly when astronauts are a part of 
the mission, the safety hazards involved must be thoroughly studied. 

Not only must precautions to prevent hazardous situations be taken, but 
also remedies must be considered if these situations occur in spite of 


The condition that the trajectories be sjnnnmetric is only a sufficient 
condition to insure repeatability. 
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the precautions. Certain unmanned experiments should be tried before a 
man is placed in the tethered body. Some effort is now being made by 
the Martin Company to detemine the feasibility of making an actual 
flight experiment to test these soft tether techniques. 


TABLE I 

Drag Acceleration Data 



Workshop 

ATM 

Coefficient of Drag 

10.15 

3.77 

Reference Area 

33.478 m^ 

33.478 m^ 

Mass 

28495 kg 

7394 kg 

m/CjjA 

83.86 kg/m^ 

58.58 kg/m^ 




Separation (km) 



Time (hr) 


FIG. 1. TYPICAL SEPARATION GROWTH BETWEEN 
TWO SATELLITES IN LOW EARTH ORBIT 
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FIG. 5. CONTOURS OF EQUAL i.,. ON THE PLANE (At.,, )f.,j 
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FIG. 6. CONTOURS OF EQUAL 





(UI) !U 



16 








17 


FIG. 8. CONTOURS OF EQUAL A ta ON THE PLANE (At^, y.,j) 
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CONTOURS OF EQUAL At, ON THE PLANE 
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CONTOURS OF EQUAL At, ON THE PLANE (At 
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FIG. 12. CONTOURS OF EQUAL TOTAL CYC 
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FIG. 13. 
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24 


FIG. 15. CONTOURS OF EQUAL VERTICAL OSCILLATIONS 
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FIG. 16. ZERO OSCILLATION CASES 
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ATION IMPULSE MADE 3 MINUTES EARLY 
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APPENDIX A 


Solution to the Linearized Differential Equation 


Equation (1) may be integrated once to yield 


X ® Dt + 2wy - j c^, (Al) 

where is a constant of integration. Substituting equation (Al) into 
(2) and rearranging yields 


y + 


P 

CO y 


■j coc^ - 2co Dt. 


(A2) 


The solution to the homogeneous part is 


yj^ * -b^ cos (jjt + a^ sin (jt, 


where bp and Cp are arbitrary constants. By inspection, the solution to 
the nonhomogeneous equation has the form 


y = Yh + Cl + Cat. 

By differentiation and substitution into equation (A2), we obtain 


yjj + co^yj^ + u)^i + to^Cat = coc^ 


2co Dt. 


Therefore, 


Cl 
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and 


Cg - 


2D 

w ’ 


2c 

y = -b cos tot + a sin wt + ~ t. 

o O JO) Cl) 


(A3) 


Substituting this into equation (Al) and integrating yields 


X = -2b sin cot - 2a cos cot + c t - -r Dt^ + d , 
o o o 2 o 


(A4) 


where do is the constant of integration. By differentiating equations 
(A3) and (A4) with respect to time, x and y are obtained as functions 
of time: 


X = -2wb^ cos cot + 2coa sin cot + c - 3Dt (AS) 

o o o 

y = ci) sin cot + coa cos cot - ~ . (A6) 

O O Cl) 


The constants of integration may be expressed in terms of the conditions 
of position and velocity at time t = 0, x^, y^, and y^: 


a 


o 



b = 
o 



c = - 3x 

o o o 


d 

o 


x + 2a^. 
o o 


32 



APPENDIX B 


The Conditions Necessary for a Repeatable Sequence 


To discover what must be done to obtain' a repeatable sequence, all 
of the equations which must be simultaneously satisfied must be written 
down. It turns out that the equations simplify somewhat if the time, 
t = 0, is placed half way into the initial trajectory; i.e., if At^ is 
the time spent in the initial trajectory between pulls, then the trajec- 
tory begins at time, t » -1/2 Ati, and ends at t = 1/2 At^. The initial 
conditions are as follows: 


'li o 


2c 




^li “ 3^7 S ■ ^o I ‘^1 " 1 


cos (lAti) 
(Bl) 

(B2) 


3 11 

X . = c + *o DAti - 2(Jb cos *;r - 2coa sin — 

11 oz o / o 2 


(B3) 


li 


- ^ - 


w 


cob sin -T ciAtj^ + toa cos -r coAtn, 


(B4) 


where Sq, b^, and c^ are defined as before. 

The final conditions just before the impulse which initiate the 
transfer are 


13 1 1 

(B5) 

2c D 1 1 

^if “ 3(? " w ■ ’^o 2 ®o 2 
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(B7) 




= c 


3 11 

THAtj, - 2(1*^ cos (oAti + 2cja^ sin ^ 


y 


If 


2d . , . 1 . 

— + cob sin T + ciH cos 

CO o 2 o 



(B8) 


At this time, an impulse is made which changes the velocity by AV^ 
in the direction of the tether, making the initial components of the 
velocity of the transfer conic 


21 If 


AVi 


^If 


(B9) 


^2i = ^If ’ 


If 


Kt + 


(BIO) 


For the transfer conic, if time is referenced to this point, after 
some time interval, Ats, it is desired that the satellite be back at the 
initial position of the first trajectory. Again, by using equations (A3) 
through (A6) , 


X 


li 


y 


li 


X 


If 


2c' 
o 

3co 


3 

+ c^Atg - "2 DCAtg)^ - 2b^ sin CkAtg + 2a^(l - cos 


- — Ato - b' cos (*Ato + a' sin oAto. 
CO o ^ o 


(BID 


(B12) 


The final velocities in the transfer conic will be 

ic - = c' - 3DAto - 2cii)' cos coAto + 2coa' sin coAto (B13) 

2f o o ^ o 

Otv 

y- = i^ + cob' sin (Atg + coa' cos cAta- (B14) 

2f CO o o 
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In the last four equations. 


a' 




= - 3x^,. 


21 


The initial velocity now must be obtained by making a second tug and 
changing the velocity by an amount AV 2 . 
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(B15) 


• ^li 

■ ^2 — 


(BI6) 


After substitution of the definitions for ao, bo> Cq, ao, b^, and Cq, 
there are 16 equations which must be satisfied by 20 variables. These 
are Xq, yo, Xq, yo, Ati, x^i, yii, x^i, y^, Xif, yif, Xif, yif, AVi, 
xai> Yaij ^^ 2 > x^f, Yaf* ^ 2 * Since there are four more variables 
than equations, four of the variables may be treated as parameters. 

Let these parameters be Xi^, y^^, yo, and Xq. If, arbitrarily, and 
Xq are set equal to zero, then bp = = 0. From equations (A3) through 

(A6) , the equations of motion then become 


Xi = 2a^ (1 


cos Ojt) 



+ X 


o 


(B17) 


Y-, = a sin t 
o 



CO 


(B18) 
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xi = 2(ja^ sin ojt - 3Dt (B19) 

2d 

yi = oSq cos tot — ^ . (B20) 


As can be seen, this makes the trajectory symmetric about the horizontal 
axis; i.e., x(t) = x(-t) and y(t) = -y(-t). Also notice that x(t) = -x(-t) 
and y(t) = y(-t). Therefore, the initial and final conditions at the 

times -1/2 At^ and 1/2 At^, respectively, are also symmetric: x^f = Xi^, 

yif = -yii, xif = -Xii, and y^f = yxi* By solving equation (Bl) for Xq, 
and equation (B2) for the variables Xq and may be eliminated from 
the system of equations, if the resulting expressions are substituted 
into the other 16 equations. Noting that bg and Cg are now zero, equa- 
tions (B3) and (b 4) become 


X . = 2oy . - T DAti 
11 11 2 ^ 


(B21) 


coAti 


+ D-^tx) cot 


- 2^ 

CO 


(B22) 


These two equations become the initial conditions which cause the motion 
to be symmetric. By using the symmetry conditions, x^f, yif, x^f, and yi^, 
can be eliminated from equations (B9) through (B12): 


X . 

21 



AV 


11 


[x^. + y^. 

11 11 



+ AVx 


[x-i + y-,]^/- 


(B23) 


(B24) 


0 = c' At 2 - 4 D(At 2 )^ - 2b' sin oiAt 2 + 2a' (1 - cos ojAt 2 ) (B25) 


2c 


, = — 2. _ -2. - b' cos GoAt 2 + a' sin cjAt 


11 w w 


2 > 


(B26) 
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where 



c* = -6wy . - 3x„. . 
o •' ll 2 l 


Using equations (B13) and (Bl4) to eliminate and j'sf equa- 

tions (B15) and (B16), 


c' - 3DAt2 - 2(jJj* cos (*iAt 2 + 2u a' sin - aV^ 


^i 


li o 


y,i “ ■ + wh' sin aiAto + coa' cos uAto - AVo 

iX Ou O O 


l“u " 


Hi * 


2 lx/2 


(B27) 

(B28) 


Now there are eight equations (B21) through (B27), with the two param- 
eters, xxi and yn, and eight variables, x^, y^i, At^, xsi, AVi, ysi, 

At 2 , and AV 2 , having eliminated the others. Equations (B25) through 
(B28) insures that the initial position and velocity components are 
reproduced. Now, Instead of making AV^ of such a magnitude that the 
initial conditions are reproducible, consider the consequences of making 
it so that the new (transfer) trajectory is symmetric about the horizontal 
axis. Because of the symmetry after some time At 2 » the body is back at 
the initial position of the first trajectory with velocity components 
symmetric to those of equations (B23) and (B24): 


2f 


-X 


2i 


X + AVx 7- 


li^ 
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y^f “ 2 ^ 2 jl/2 

" 11 -'ll 


After the second pull on the tether, the velocity components are 


X , X 

X = X, . + AV, — AVo 


. , ’’ll ... 

y = y, . + AVn AVo 


Notice that if the second pull is made so that AV 2 =* AV^, then 


X = X , 

new 11 


yn^ry “ ^ ^ 
new 11 


i.e., the initial velocity has been reobtained. Therefore, the equations 
(B25) through (B28) insuring that the initial conditions are reobtained 
can be replaced by making AV^ so that the trajectory is symmetric, wait- 
ing some time Atg, and making AV 2 = AV^. To insure symmetry, AVi must 
be made so that Xq^, and y 2 i in equations (B23) and (B24) must satisfy 
equations similar to equations (B21) and (B22) after the time At 2 : 


Xgi - ” 2 15At2 - " 2 ^^^2 


(B29) 


y2i = — r 


= ((^ 3 ^. + DAt 2 > cot ~Y 


(B30) 
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Notice that the four equations (B25) through (B28) have been replaced 
by equations (B29) and (B30) and the condition AV 2 = AV^. This indicates 
that there was some redundancy in the system of equations (B21) through 
(B28), and there is one more variable which may be arbitrarily treated 
as a parameter. Let it be At^. Before parameterization, more variables 
and equations can be eliminated by replacing AVs by AV^ and by eliminat- 
ing xsi and ysi between equations (B23), (B24), (B29), and (b 30). The 
system of equations then becomes 


^li " ^‘^li ■ 2 


wAtj. 

y^i - + DAti) cot ^ - 2 - 


[x^i + 


+ AVj 


+ DAtg) cot — y- - 2 ^ 


[xfi t 


Eliminating and yields 


AVn 


^li 




= 2 D(Ati + Ata) 


(B31) 


AVi 


11 




icAti 

(wy^i + DAtg) cot — ^ + (wy^^ ■ DAti) cot — — 


(B32) 
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Eliminating AV^ yields 


j D(Ati + At 2 ) - — = + DAta) cot — j— + (wy^^ - DAtj.) cot -y 


(B33) 

Thus, all of the pertinent variables can be expressed in terms of the 
parameters xi^, yii, and At^. Equations (B21) and (B22) give and 
yii in terms of y^i and At^. Equation (B33) can be numerically solved 
for Ats. After this is done, AV^ can be found from equation (B31), 
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